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Abstract 

Tunneling in quantum field theory is well understood in the case of a single scalar field. 
However, in theories with spontaneous symmetry breaking, one has to take into account the 
additional zero modes which appear due to the Goldstone phenomenon. We calculate the 
false vacuum decay rate in this case. The result may differ significantly from the tunneling 
rates in the absence of symmetry breaking. 
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The effects of tunneling between different vacua in quantum field theory play an important 
role in high-energy physics, cosmology and condensed matter physics. The rate of the false 
vacuum decay was calculated in Refs. |], 0, [| for a potential which depends on a single 
scalar field. However, most theories beyond the Standard Model, e.g., the MSSM, grand 
unified theories, and strings, employ scalar potentials which depend on several scalar fields, 
and may have degenerate minima related by some internal symmetry. We will see that in 
this case the results of Refs. [|T], |2], ||] do not apply because of the complications that arise 
due to the additional zero-modes associated with internal symmetries. The purpose of this 
letter is to generalize the semiclassical calculation of Refs. [§, |[] to such cases. 

If the potential U((f>) depends on a single scalar field <p(x), then the corresponding path 
integral, 



[d<j>] e 



-S[<t>] 



is dominated by the field configuration (f>(x) called the "bounce" and can be evaluated using 
the saddle point method || |3|. Then the transition rate per unit volume, in the semiclassical 
limit 01 is 



r/v = l(s[0])V s[<s/?i 



det'[-# + U"{4>)] 



-1/2 



x(l + 0(h)) (2) 



det[-<92 + [/"(0)] 

where S[cj)} is the Euclidean action of the bounce, and det' stands for the determinant with 
the four zero eigenvalues omitted. These four zero modes are related to the translational 
invariance in the four- dimensional Euclidean space. Integration over each of them yields a 



factor of yS[(j)} in equation (H). 

Let us now consider a quantum field theory with a scalar potential U(<pi, ...,4> n ) which 
has a local minimum at <pi = <p2 = ■■■ = 4>n = 0, Z7(0, 0) = as well as at least one 
additional (local, or global) minimum at 4>i — 4>ti i = 1)2, ...,n; U{(f>\, ...,(j)n) < 0- 

The bounce, (f)(x) = (<pi(x), <p n {x)), is the stationary point of the action, so that 
SS[4>] = 0. It is found as an 0(4)-symmetric || solution 0(r), r = ^Jx~^xP, of the correspond- 
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Figure 1: The potential U(<f)i, 02) with several degenerate minima related by a Z m symmetry 
ing Euler-Lagrange equations^: 

^i(r) = 4r U (^-^n) (3) 
with the following boundary conditions which ensure the fmiteness of S[(f>]: 

' (d/dr)Ur)\r=o = Q 

< _ (4) 
^ 0i(oo) = 

If the potential U (0) has only one local minimum characterized by a negative value of the 
potential, equation (H) gives the correct answer for the rate of the false vacuum decay. Let us 
now consider the potential U(<f>i, 02) which has m degenerate local minima related by the Z m 

1 In most cases, one cannot solve this system of non-linear equations analytically. An effective method for 
finding the bounce was proposed in Ref . || . 
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symmetry as shown in Fig. 1. The equation (0) will have m independent solutions satisfying 
the boundary conditions each bounce pointing in the direction of a given minimum in 
group space. Tunneling to different vacua are independent events, and therefore the total 
probability of the false vacuum decay is the sum of the contributions from each individual 
bounce. The naive application of the formula (0) might lead one to conclude that as the 
number of local minima m increases, the lifetime of the false vacuum would approach zero. 
If one then considers the case of a broken U(l) symmetry (Fig. 2) as the limit m — > oo, one 
might argue that the false vacuum would decay infinitely fast because there are infinitely 
many channels, each having a finite probability determined by formula (0). This conclusion 
is wrong because, in the case of a continuous symmetry, the Goldstone modes give rise to 
some additional zero modes in the determinant in equation (^), so that the right-hand side 
of (0) is, in fact, divergent. 

Clearly, formula @ cannot be used to describe the tunneling probability in theories with 
spontaneously broken continuous symmetries. One must not include the Goldstone zero 
modes in the determinant, but deal with them instead in a similar manner to the way one 
usually treats the zero eigenvalues associated with translational invariance. The latter is the 
basic idea of our calculation. 

We begin by considering the potential U (0i, 02) = ^(0i + 02) which has a U{1) symmetry 
(Fig. 2). We suppose that in addition to the local minimum at = 0, U((f>) has a U{1) orbit 
of degenerate global minima at some non-zero value of 0. It is convenient to redefine the 
fields in terms of the radial and the angular components: 



Then U((f)) = U(r)), (d/d6)U = 0, and every solution of the equation (^) can be written 
in the form 




(5) 




,6 



(6) 
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Figure 2: The potential £/(0i,0 2 ) = U{<p\ + 0|) with a U(l) orbit of degenerate minima. 



where 9 is a constant independent of x and fj(x) is the usual bounce, the same as in the 
one-component case. 

The idea of the method, used in Refs. |2|, |J to calculate the transition probability, was to 
evaluate the path integral using the saddle point method. The bounce <f)(x) is a stationary 
point of the action (<5>S[0] = 0), but it is not a minimum because the second variation 
operator 

ss _ A | d 2 U(<p) 
5(f)i5(f)j lj d(f>id(f)j 

has one negative eigenvalue. In addition, this operator has four zero eigenvalues due to trans- 
lational invariance in the coordinate space. The corresponding normalized eigenf unctions are 
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<p, = {s$\y l / 2 d^{x), /i = i,...,4 (8) 

In the vicinity of the bounce, one can expand the field 4>(x) in the (complete orthonormal) 
basis of the eigenvectors of the operator (|7|) : 

fe (x) = fc (x) + $>M } (9) 

i 

Then the path integral measure, [dcf>k], is defined as 

[d<l> k ) = Y[(2ich)~ h de < £ ) (10) 

i 

The path integration with respect to the Fourier coefficients cjj. , which correspond to 
positive modes, is Gaussian and results in the determinant of the type (0) where the zero 
modes are not included. The integration over zero modes, however, is less straightforward. 
We would like to treat the Goldstone modes on the same footing as the Poincare zero modes 
(§). The first step is to find the orthonormalized basis in the subspace of eigenvectors that 
correspond to zero eigenvalues. By explicit differentiation of equation (|3|) with respect to 9 
we find that, as long as U((f>) = U(rj) is independent of 9, 

(j )e = (| fj 2 (x)d*x) 8e4>{x) (11) 

is an eigenvector of (|7]) corresponding to a zero eigenvalue. The normalization comes from 
/ |9e0(a;)| 2 (i 4 x = / fj 2 (x)d 4 x. It is easy to see that the eigenvector (|HD is orthogonal to 
each of the eigenvectors @. Therefore, the set of functions (p|) and ( II ) combined form a 



complete orthonormal basis in the subspace of zero modes of the operator (0). 

Let ce be the Fourier coefficient in (|9|) corresponding to the Goldstone mode. Then the 
change in <p corresponding to some small change in eg is d(f) = (pgdcg. On the other hand, 
d<p = dg<f) d9. Comparing these relations to each other, and using the normalization from 
equation (|TT|), °ne gets 
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(2Txh)- l l 2 dc e = ( / ff{x)d A x/2TiTi) 1 l 2 de (12) 
Therefore, integration with respect to cq in the path integral ([!]) is equivalent to the 



integration over the angle 6 in group space, up to the factor of J J ff. The field configuration 
that comprises n widely separated bounces is also a saddle point of the path integral with the 
action nS[<j)}. We have to sum over all such configurations and integrate over the positions 
of the centers of the bounces. This is done in full analogy to the corresponding calculation 

in 1,1- 

The result is 



r/V = 2tt 



r] 2 (x)d 4 x 



1/2 fS[fj]' 



-s[fj]/h 



2nh 



det"[-d 2 + U"(f))\ 



det[-d 2 + U"(0)\ 



-1/2 



x(l + 0(h)) (13) 



where det" stands for the determinant from which all zero eigenvalues, including the Gold- 
stone zero mode, are excluded. 



The result (13) can be easily generalized to the case of a larger internal symmetry group. 
Suppose the true vacuum is characterized by the vev that breaks some symmetry group G 
down to G', so that the scalar potential has iV = rank(G/G') Goldstone modes. The false 
vacuum decay rate is then 



r/v = C{G/G') 



r\ [x)d x 



N / 2 fS[f}\ 



-s[ff\/h 



2nh 



det"[-dl + U"m 



det[-8g + U"(0)] 



-1/2 



x(l+0(h)) (14) 



Here C{G/G') is a dimensionless factor that depends on the quotient group G/G' . The 
pre-exponential factor in square brackets has dimension N in length units, and the determi- 
nant in the numerator, with 4 + N zero modes omitted, has dimension 2(4 + N) in excess 
of that of the denominator. The power (—1/2) reduces the latter to —(4 + N), so that the 
total dimension of the right-hand side is minus four, as it should be. 

Some comments are in order. We note that the result ( |14|) can differ significantly from 
that obtained by using the formula (E3) and ignoring the symmetries of the potential. For 



6 



example, let us consider the case of G = SU(5) — > G' = SU(3) x SU{2) x U{1). There are 
N = 12 independent Goldstone modes in this case. In the thin wall approximation, one can 
express / ff in terms of S[fj\ times some dimensionful quantity. Then it can be shown that 
the answer obtained by applying formula (|2|) and ignoring the Goldstone modes would be off 
by a factor of (>S[0]) 6 . For S[4>} ~ 10 — 100, the difference would be as great as six to twelve 
orders of magnitude. 

We also note that the increase in the number of degenerate minima m (Fig. 1), does not 
necessarily imply an increase in V. In fact, to accommodate m degenerate minima, m — > oo, 
separated by the fixed height barriers, one has to allow the curvature (in the ^-direction) 
at the bottom of each minimum to grow as ~ m 2 . This corresponds to ~ (1/m) behavior 
of the determinant factor in equations ( |T3"D or (0). Now the total decay rate is m times 
the contribution of each minimum, V ~ r o m(l/m) — > const. This is the resolution of the 
paradox posed in the beginning. 

It is straightforward to generalize formula ( [Hp to the case of quantum field theory at 
finite temperature 0. 

We have assumed that G' C G, which is the case for many applications in particle physics. 
However, in some cases the true vacuum may, at some temperature, have a larger symmetry 
than the false vacuum (see, e. g., Ref. [0]), so that G C G'. Finally, the true vacuum may be 
characterized by the symmetry group, G', completely unrelated to that of the false vacuum, 
G. If this is the case, then in general one has to modify the expression (0) to take into 
account the particular features of the potential. 

To summarize, we have calculated the tunneling probability in quantum field theory with 
spontaneously broken local or global symmetries. 
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comments. We are also grateful to G. Segre and P. Steinhardt for useful conversations. This 
work was supported by the U. S. Department of Energy Contract No. DE-AC02-76-ERO- 
3071. 
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